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PREFACE 


A new technique has been developed for the weighting of data from 
satellite tracking systems in order to obtain an optimum least squares 
solution and an error calibration for the solution parameters. Data 
sets from optical, electronic, and laser systems on 17 satellites in 
GEM-T1 (Goddard Earth Model, 36x36 spherical harmonic field) have been 
employed toward application of this technique for gravity field 

parameters. Also GEM-T2 (31 satellites) was recently computed as a 
direct application of the method and is summarized here. The method 
employs subset solutions of the data associated with the complete 

solution and uses an algorithm to adjust the data weights by requiring 
the differences of parameters between solutions to agree with their 
error estimates. With the adjusted weights the process provides for an 
automatic calibration of the error estimates for the solution 

parameters. The data weights derived are generally much smaller than 
corresponding weights obtained from nominal values of observation 
accuracy or residuals. Independent tests show significant improvement 
for solutions with optimal weighting as compared to the nominal 

weighting. The technique is general and may be applied to orbit 
parameters, station coordinates, or other parameters than the gravity 
model . 
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INTRODUCTION 


The method of data weighting has been an outgrowth of a 
calibration process for the error estimation of gravitational models 
that have been derived from satellite data, Lerch et al. (1985 and 1988 ) 
and Wagner and Lerch (1978). The principle of the new technique is to 
estimate the weighting of the data so as to produce realistic error 
estimates of the solution parameters from subset solutions of least 
squares normal equations. Application has generally been with use of a 
large set of satellites with inhomogeneous data from tracking systems of 
laser, electronic, and camera (optical) data. The gravity model of 
GEM-T1 (Marsh et al., 1988) using some 17 satellites has been tested 
with the new technique and the GEM-T2 (Marsh et al . , 1989) solution with 
some 31 satellites has been derived with the new process of optimum 
weighting of the satellite data sets. 

The accuracy estimation of the gravity model is particularly 
important for the TOPEX Project (1992 launch) for ocean application of 
its altimetry. It requires that the radial orbit error be accurate to 
better than 10 cm due to the uncertainty of the gravity field. Hence 
the estimation process for the errors, which are based upon the weights 
assigned to the data, must be reliable. The accuracy of the solutions, 
particularly the low degree field, is also important for the Lageos 
orbit. Accuracy is needed for the estimation of baseline motion of 
laser tracking sites at the centimeter-per- year level as part of the 
NASA Crustal Dynamics Project. 
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2 . 


OBSERVATION WEIGHTING AND DATA CHARACTERISTICS 


Observations obtained from geodetic satellite tracking systems 
generally have precision levels, particularly laser systems, much better 
than the observation residuals obtained from satellite orbital arcs in 
post-fit analyses using the best models. This is true even though the 
orbital models employed were derived from the same satellite data and 
with the same arc lengths of several days. The problem here is that 
there are unmodeled systematic errors (biases) which need to be 
accounted for in the weighting system of the solution (Brown, 1988). 

In Figure 1 an example of the characteristics of the residuals is 
shown for a pass of data from a typical laser tracking site. The 
precision error (a ) of the laser data is generally small (centimeter 
level) as compared to the rms (o fc ) of the residuals for a satellite data 
set t. Values of are given in Tables 1A and IB (GEM-T1 and T2 data 
sets) for different satellite data types, and for laser systems, o fc varies 
from 10 cm for Lageos orbits to over 50 cm on GE0S-1 orbital data in 

1978. 


Note in Figure 1 that the residuals of a tracking pass with noise 
removed fit very closely to a straight line as a function of a bias 
offset (b ) and a timing error. The bias offset is the dominant part of 
the residuals. If the residuals were random with rms equal to the 
weight per observation point should be 



but with a constant offset (b Q ), say for N=50 points in a pass, the 
weight should be degraded by 

w^ = 1/No^ = .02/o^ . 
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The latter case is characteristic of our situation particularly for the 
high-precision laser data. The bias effects per pass tend to fluctuate 
randomly from pass to pass. 


In general, for a given satellite data type t we have 


w 


t 


v° 


2 

t 


where cr is the rms of residuals for the satellite data set and f.. is a 
t u 

downweighting factor to account for the bias effects and the correlated 
effects of the residuals, particularly within the pass. The weighting 
technique will obtain w t directly. Note from Table 1A (and IB), o t as 
well as Oj. is given for each data type where 



hence 


f t = H /o t )2 

which is approximately a constant 
f t - .01 

for the satellites with the laser data. In Table IB for the Starlette 
( 1 86 ) and AJISAI laser data.f * .002 instead of .01 since the data 
rates were 5 times faster (1 per second as compared to 1 per 5 seconds), 
and hence f t counteracted the excess data. Note also for the optical 
systems, the factor f t is generally much larger than with the laser 
systems. This may be expected since the number of points per pass for 
optical data are fewer and the ratio of noise to bias is much larger 
than with the laser data. 
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LEAST SQUARES MINIMIZATION 


The method of solution is a modified least squares process which 
minimizes the sum (Q) of signal and noise as follows: 


Q 


1 


m 


i ,m 


i I y't 

t obs a t 
i 


(1) 


where the signal is given by 


spherical harmonics comprising the solution 
U i,m’ coefficients; and 

is rms of the coefficients of degree a (a priori 
rule) and is introduced to permit larger solu- 
1 10" tions to degree and order 36x36. This law, 

°sT - x 2 based upon Kaula's rule, has been obtained inde- 

- 1 pendently from studies of the spectra of the 

Earth's gravity field and is used here to repre- 
sent the observed power within the geopotential. 


and the noise by 

observation residual (observed-computed) 
r : for the ith observation of satellite 

il ' tracking data set (type) t; and 

rms of observation residuals (generally 
o fc : significantly greater than a priori 

data precision) 

downweighting factor to compensate for unmodeled 
*t * error effects for each data type t (ideally f=1 for 
pure noise). 


The optimum weighting method estimates the combined weight 
directly, namely 



( 2 ) 


When minimizing Q above using the least squares method, the normal 
matrix equation and error covariance is obtained as follows: 

are the normal equations, where x is the 
solution, R is the vector of residuals, and 

is the approximate form for the variance- 
covariance error matrix which must be 
calibrated by adjusting the weighting. 

The process of minimizing both signal (Kaula constraints) plus 
noise in (1) is also known as collocation by Moritz (1978). With the 
normal least squares approach (noise-only minimization) there is a 
problem of separability due to the strong correlation between many of 
the high degree coefficients. The absence of collocation (GEM-T1 without 
the Kaula constraint) results in excessively large power in the 
adjustment of the potential coefficients. Figure 2 illustrates the 
instability of the least-squares solution when collocation is not 
used. A satellite-derived gravity solution has been solved without 
collocation which is evaluated using a global set of independent gravity 
anomalies. An unrestricted high degree field performs poorly due to 
excessively large adjustment in the coefficients which is normally 
circumvented in the standard least-squares method by solving for a 
smaller sized field. Unfortunately, by restricting the size of the 
field, one also is requiring the higher degree terms above the field 
limits to be constrained absolutely to zero. Figure 2 also shows the 
disadvantage of this approach where the smaller sized field (PGS-3067) 
contains aliasing in its coefficients and does not perform well. (The 
abbreviation PGS stands for Preliminary Gravity Solution.) The aliasing 
signal sensed in the data above the field limits is absorbed into the 
adjustment of the lower degree coefficients. The best approach is seen 
with the least squares collocation (or constrained) solution, GEM-T1, 
with a complete solution of a 36x36 field in harmonics.. 


Nx = R 

V = N~ 
zz 


4. 


LEAST SQUARES NORMAL EQUATIONS 


In matrix form the observation equation is given by, assuming 
linearity, 


0 - C = r = r - Ax 
o 

(3) 

x = X - X R 


where 


r = 0 - C — residual, observation (0) minus computed value (C) 
from solution 


x = X - X R adjustment of solution (X) from reference value 

( Xpj ) (for error analysis X R = X(true)) 

A — matrix of partials evaluated at X = X R 


r 


o 


— residuals based upon a priori value X R . 


For the gravity field, the linearity of perturbations may be seen for 
the spectrum of harmonics in Kaula (1966). The weighted normal 
equations are where W is a diagonal weight matrix (Lawson and Hanson, 
1974) — 


A T Wr = 0 
then from (3) 

A^WAx = A^Wr . (4) 

o 

For error analysis it is convenient to let the reference value 
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X R = X( true) 


then from (3), x is the error in the solution X, namely 
x = X - X(true). 

Hence (4) becomes 

A T WAx = A T We (5) 


where 


C ' r 0 ■ 0 - C B 
= 0 - C(true) 


represents the errors due to all unmodeled systematic effects including 
random noise but excluding errors in the adjusted parameters. Instead, 
these are the errors in x given by the solution to (5). Our solutions 
will be represented by the form (5) as we are interested in the 
difference between two solutions, x and x, namely 


x - x = [X - X(true) ] - [X - X(true)] 
= X - X . 


( 6 ) 


The normal matrices for (5) are written compactly as 
Nx = R 


where 
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( 7 ) 


N = A WA 


R = A We . 


The normal matrices for each data subset t will be given as 


"t “t = "tW 


( 8 ) 


“A * “A e t 


where t=0 is a special case which corresponds to the signal constraints 
where the weight is fixed. 
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5 . 


METHOD OF ESTIMATION OF WEIGHTS 


The technique for estimating w fc for each data set t is based upon 
a Complete solution (S) with all the data and a subset solution (S^) 
where data set t is removed. Let the normal equations for the complete 
solution x and the subset solution be given as in (7), namely 

Nx = R (g) 

N x t = R 

where from (7) 


N = l w N 
J*t J J 


R 


l 

j*t 


W J R J 


N = N + w^N^ 


R = R + w t R fc . 


( 10 ) 


The covariance (variance-covariance) matrices (V) for the errors x 
and x^. are obtained as 

V(x) = N ^ = E(xx^) 

(11) 

V(x fc ) = N _1 = E(x fc xJ) . 

As in (6), 


x fc -x = [X t - X(true)] - [X - X(true)] 


The covariance of the difference between the solutions is 


( 12 ) 
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(13) 


T 

V(x fc -x) = E(x t -x)(x t -x) 

= V(x fc )-2 E(x fc x T )-V(x) 

= V( x fc ) - V(x) 

where as shown below 

E(x t x T ) = V(x) . 

From (9), 

E(x t x T ) = N" 1 E(R R T ) N' 1 
From ( 10) and (11), 

E(R R T ) = E[R(R + w t R t ) T ] 

= E(R R T ) = H V(x fc ) N < l6 > 

= N 


since 

E(R R^) = 0 . 

The latter result is true as from (10) the data set t is excluded from 
the subset solution, making R and R fc independent. Hence (14) results 
by substituting (16) into (15) and using (11). 
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5.1 WEIGHTING ALGORITHM 


Using just the gravity parameters in (X t -X) the weighting 
algorithm is given by the calibration factor k fc obtained from 


(X t -X) T (X t -X) = (x t -x) T (x t -x) = k fc TR V(x t -x) 


(17) 


where TR denotes the trace of the matrix and where from (9) through (13) 


x t -x = N _1 R - N -1 R = X fc -X 


N = l w N 
J*t J J 


( 18 ) 

(19) 


N = N + w N 


l w j R j 

+ f- J J 




( 20 ) 


R = R + w t R t 


V(x t -x) = N 1 - N = V(x ) - V(x) 


( 21 ) 


Since k t scales the error variances it will be inversely 
proportional for scaling the weight w^. to obtain the adjusted weight 
w£, namely 


w t = w t 7 k t • 


( 22 ) 


This latter result will be derived more directly below. By iterating on 
the solutions x^ for each data set t and the complete solution x for all 
data sets until 


k t = 1 
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for each t, the weights by (22) will then converge and the error 
estimates will automatically be calibrated from (17). 

Results are given below to show how the weights and associated 
calibration factors converge. Because of the extensive computations for 
a large number of data sets a reasonable set of a priori values for the 
weights should be available for their refinement in the optimization 
process. 

The gravity parameters of spherical harmonic coefficients are 
calibrated as a set by (17). Calibrations (k t ) are also given by 
subsets of spectral components from the harmonics of degree l and order 
m. For all satellite data sets t (Lerch et al., 1988) relatively little 
variation is seen in the spectral calibrations. 


5.1.1 Weighting Adjustment 

The relation (18) for the weighting adjustment 

w t = V k t 

is derived from use of (17) through (21). It is assumed that the data 
set t does not significantly change the solutions x and x^ beyond first 
order effects as follows: 

V(x^-x) = N ^ - N ^ = N ** - (N+w^N^) (23) 

= N" 1 - (I+w t N t ) _1 N _1 
* w JT 1 N fc N _1 . 

To the same approximation 
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V x 3 w t N_1 R t 

E(x t -x)(x t -x) T - w* N' 1 E(R t R^) N _1 . (24) 

From (8), 

E(R t R t T ) = aJ E(e t eJ) A t (25) 


where o t accounts for the unmodeled systematic effects in e t and the 
corresponding weighting effect is given as 


t " "2 ~ t t ♦ 

°t 

Using (23) and (25), then (24) becomes 

T W t 

E(x t -x)(x t -x) = ~ V(x fc -x) 
From (26) and (17) 

k t = W t /W t 

which gives the result (22). 


( 26 ) 
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6 . 


TESTS AND RESULTS FOR OPTIMUM WEIGHTING TECHNIQUE 


Sample tests of the weighting algorithm (22) were made using 
GEM-T1 plus additional data sets for several satellite data types of 
laser, optical, and electronic data. Results are given in Table 2 which 
show that the algorithm nearly converges in one step from the a priori 
starting weights. Plots of w^ vs k^. from (17) show a strong linear 
relationship from the origin (w^. = k^ = 0). Hence 

vT w 
k' = k 

and by setting k' = 1 for calibration the adjusted weight w' should 
nearly converge from (22). 

The above tests were made in preparing the weights for additional 
data sets to GEM-T1 that were combined for the GEM-T2 model. The 
convergence of these weights for GEM-T2 is shown in Table 3. In 
addition to the optimum weights the technique provides an automatic 
calibration of the error estimates based upon the satellite data types t 
since each of the k t from (17) is required to converge to 1. 

The data weights in GEM-T1 were derived primarily by requiring the 
weight for each data type t to give the best overall agreement with 
independent mean gravity anomalies (Rapp, 1986) and with the satellite 
observation residuals on selected test arcs. The calibration factors 

(k ^2) for several of the major data types (Lerch et al . , 1988) are 
t 

given in Table 4 which show that the weights converge (k fc = 1) except 
for the Lageos laser data. However, several additional tests were made 
in Table 4 for the calibration factor using independent data from Seasat 
altimetry (Rapp, 1986) and surface gravity data (Pavlis, 1988). All of 
the latter tests show good calibration of the error estimates, 
indicating optimum weighting was closely achieved. The last test 
deliberately increased the weighting for a subset of laser data by a 
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factor of 10 giving a value k t =(2.75) 2 . From (22) the adjusted weight 
should be reduced by a factor of 1/k fc which would nearly recover the 
original weight in one step of the iteration process. The gravity model 
with the increased weight naturally gives smaller error estimates but it 
also gave significantly worse agreement with independent surface gravity 

anomalies. 


15 



7. 


SUMMARY 


The optimum weighting technique was shown to be important in the 
weighting of satellite data, particularly precise laser data where 
unmodeled systematic effects require a significant downweighting factor 
as shown in Table 1. The method of weighting was shown in Section 5.2 
to provide realistic error estimates for GEM-T1 and-T2. These models 
were calibrated using subset solutions based not only on data employed 
in their solutions but also upon independent data from altimetry and 
gravimetry. Because of the important application of the gravity model 
to ocean altimetry in the Topex Project, the gravity model errors were 
projected on the radial component of the TOPEX orbit and the result gave 
10 cm for GEM-T2 which nearly meets the goal of the gravity model. 

It was also shown in Section 5.2 that the model with the increased 
weight on the data over the optimum weighting gave much poorer agreement 
with independent surface gravity anomalies. The optimum weighting 
technique based upon the mathematical formulae is general and may be 
applied to other than gravitational parameters such as station 
coordinates and, in particular , orbit parameters where knowledge of 
accuracy estimation and refined solutions are needed. 
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TABLE 1A 

SATELLITE DATA IN GEM-T1 


SATELLITE 


LAGEOS 
STARLETTE 
GEOS-3 
PEOLE 
BE-C 

6 GEOS-1 

7 GEOS-2 

8 Dl-C 

9 Dl-D 

10 SEAS AT 

11 OSCAR- 14 

12 AHNA-1B 

13 BE-B 

14 COURIER-1B 

15 TELSTAR- 1 

16 VANGUARD -2Rb[ 

17 VANGUARD-2 


SEMI MAJOR 
AXIS (km.) 


12273. 

7331. 

7226 . 

7006. 

7507. 

8075. 

7711. 

7341. 

7622. 

7170. 

7440. 

7501. 

7354. 

7469. 

9669. 

8496. 

8298. 


ECC 

INCL 

DEG 

DATA 

TYPE 

9 OF 
ARCS 

* OF 
OBS 

RMS 

RESID. 

°t 

SIGMA* 

WEIGHTS 

K 

.0038 

109.85 

LASER 

57 

144527 

10cm. 

112cm. 

.0204 

49.80 

LASER 

46 

57356 

20cm. 

224cm. 

.0008 

114.98 

LASER 

36 

42407 

70cm. 

816cm. 

.0164 

15.01 

LASER 

6 

4113 

90cm. 

816cm. 

.0257 

41.19 

LASER 

39 

64240 

50cm. 

577cm. 



CAMERA 

50 

7501 

2 arcsec 

5.6 arcsec 

.0719 

59.39 

LASER 

48 

71287 

70cm. 

667cm. 



CAMERA 

43 

60750 

1 arcsec 

8.9 arcsec 

.0330 

105.79 

LASER 

28 

26613 

80cm. 

816cm, 



CAMERA 

46 

61403 

1 arcsec 

8.9 arcsec 

.0532 

39.97 

LASER 

4 

7455 

150cm. 

816cm. 



CAMERA 

10 

2712 

2 arcsec 

7.3 arcsec 

.0848 

39.46 

LASER 

6 

11487 

100cm. 

816cm. 



CAMERA 

9 

6111 

2 arcsec 

8,9 arcsec 

,0021 

108.02 

LASER 

14 

14923 

70cm. 

707cm. 



DOPPLER 

14 

138042 

.5cm/sec 

7cm/sec 

.0029 

89.27 

DOPPLER 

13 

63098 

Icm/sec 

8cm/sec 

.0082 

50.12 

CAMERA 

30 

4463 

2 arcsec 

4.5 arcsec 

.0135 

79.69 

CAMERA 

20 

1739 

2 arcsec 

4.5 arcsec 

.0161 

28.31 

CAMERA 

10 

2476 

2 arcsec 

4.5 arcsec 

.2429 

44.79 

CAMERA 

30 

3962 

2 arcsec 

4,5 arcsec 

.1832 

32.92 

CAMERA 

10 

686 

2 arcsec 

4,5 arcsec 

.1641 

32.89 

CAMERA 

10 

1299 

2 arcsec 

4.5 arcsec 


SIGMA (<t)=(v) 
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TABLE IB 

NEW SATELLITE DATA IN GEM-T2 IN ADDITION TO GEM-T1 


SATELLITE 

SEMI MAJOR 
AXIS (km.) 

ECC 

INCL 

DEG 

DATA 

TYPE 

#OF 

ARCS 

#OF 

DBS. 

134093 

RMS 

RESID 

5 t 

SIGMA* 

WEIGHTS 

LAGEOS 
•84. '85, '86. '87 

12273 

.0038 

109.85 

LASER 

29 

10cm. 

112cm. 

STARLETTE 

'83/84 

STARLETTE 

•86 

7331 

.024 

49.80 

LASER 

LASER 

38 

73 

40041 

411102 

20cm. 

20cm. 

224cm. 

500cm 

AJISAI 

1500 

.0006 

50.0 

LASER 

36 

156021 

1 6cm. 

316cm. 

G EOS-1 '80 

8075 

.0719 

59.39 

LASER 

30 

54129 

32cm. 

258cm. 

G EOS-3 '80 

7226 

.0008 

114.98 

LASER 

50 

54526 

25cm, 

224cm. 

G EOS-3 
G EOS-3 :ATS 
’75/76 
GEOS-3:ATS 
*77/78/79 




LASER 

SST 

SST 

26 

9 

17 

17027 

19074 

8326 

70cm. 

.4cm/$ec 

.2cm/sec 

816cm. 
7.1 cm/sec 

3.2cm/sec 

NOVA 

1170 

.0011 

89.96 

DOPPLER 

16 

73238 

,4cm/sec 

2.6cm/sec 

LANDSAT-1 

900 

.0012 

99.12 

DOPPLER 

10 

26426 

1 .5cm/sec 

1 0.5cm/sec 

GEOSAT 

800 

.0008 

108.0 

DOPPLER 

13 

549141 

1 .3cm/sec 

4.5cm/sec 

OVI-2 

8317 

.0184 

144.27 

CAMERA 

4 

973 

2 arcsec 

5.8 arcsec 

ECHO-1 RB 

7966 

.0118 

47.21 

CAMERA 

32 

4482 

2 arcsec 

8.2 arcsec 

SECOR-5 

8151 

.0793 

69.22 

CAMERA 

13 

726 

2 arcsec 

5.8 arcsec 

INJUN-1 

7316 

.0079 

66.82 

CAMERA 

44 

3310 

2 arcsec 

8.2 arcsec 

TRANSIT-4A 

7322 

.0076 

66.82 

CAMERA 

50 

3832 

2 arcsec 

8.2 arcsec 

5BN-2 

7462 

.0058 

89.95 

CAMERA 

17 

820 

2 arcsec 

8.2 arcsec 

OGO-2 

7341 

.0752 

87.37 

CAMERA 

16 

1207 

2 arcsec 

8.2 arcsec 

OSCAR-7 

7411 

.0224 

89.70 

CAMERA 

4 

1862 

2 arcsec 

5.8 arcsec 

MIDAS-4 

9995 

.0112 

95.83 

CAMERA 

50 

_ 

31779 

2 arcsec 

8.2 arcsec 
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TABLE 2 

TEST FOR OPTIMUM WEIGHTING TECHNIQUE WITH GEM-T1 AS 
SUBSET SOLUTION (TWO ITERATIONS) 


w' 


w 


t 


t " k 


t 


GEM-T1 + 


1980 GEOS- 1 LASER 
(30 ARCS) 


STARLETTE LASER 
(73 1986 ARCS) 


NOVA DOPPLER 
(16 ARCS) 


9 NEW OPTICAL SATS. 
(230 ARCS) 


LANDSAT S-BAND 
(10 ARCS) 





.49 

.05 

.10 

.88 

.10 

.11 

.46 

.020 

.043 

.78 

.043 

.055 

1 .60 

.1 

.062 

1.02 

.062 

.061 

3.2 

.2 

.063 

.97 

.063 

.065 

.60 

.0025 

.0042 

.98 

.0042 

.0043 
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TABLE 3 

DATA WEIGHTS AND CALIBRATION OF GEM-T2 


SUBSET 

SOLUTION 

DATASET 

PGS3429 

CALIBRATION 

FACTORS 

PGS3429 

WEIGHTS 

PGS3454 

WEIGHTS 

PGS3454 

CALIBRATION 

FACTQBS 

PGS3480 

WEIGHTS 

PGS3480 

CALIBRATION 

FACTORS 

GEM-T2 

WEIGHTS 

GEM-T2 (2) 
CALIBRATION 
E5CIQBS 

AJISAI 

1.28 

A 

A (,) 

1.21 


1.29 

J. 

.79 

LAGEOS 

1.29 

.8 

.8 

1.00 

.8 

1.11 

.8 

.87 

STARLETTE 

1.04 

.2,. 2, .04 

.2,. 2, .04 

1.01 

.2, .2, .04 

.96 

.2, .2, .04 

.96 

4-LASER* 

1.02 

.015 

.015 

1.00 

.015 

.96 

.015 

1.01 

GEOSAT 

.59 

.01 

■01-5. 

.66 

.035, 

.75 

JL5. 

.81 

GEOS-3ATS 
LASER, SST 

.68 

.015, .1,. 02 

■015..05..02 

.73 

.015.J_.-02 

.66 

.015,.1 ,.02 

.66 (3) 

NO/A 

.82 

.07 

■ 075 

.83 

A. 

.83 


.90 

LANDSAT 

.90 

.0075 

.0075 

.90 

.009 

.92 

.009 

.92 

1986 GEOS-3 
LASER 

.86 

.1 

JL5. 

.91 

xi 

.97 

.2 

.96 

1980 GEOS-1 
LASER 

.87 

.1 


.97 

.15 

.99 

.15 

1.05 

OPTICAL* 

.95 

.05, .06 

.05, .06 

.95 

.05. .06 

.94 

.05, .06 

.92 

SEASAT 


.02 

.02 

1.02 

.02 

.97 

.02 

.94 

OSCAR 


.015 

.015 

1.47 

JLQJL 

.95 

.007 

1.13 

3-LASER* 


.015 

.015 

.82 

.015 

.83 

JL2. 

.87 


1 . UNDERLINED WEIGHTS ARE THE ADJUSTED ONES IN THE ITERATED SOLUTIONS 

2. CALIBRATION FACTORS ARE CONSERVATIVE BUT SUFFICIENTLY CONVERGED 

3. ATS SST WEIGHT DEUBERATELY UNDERWEIGHTED BASED UPON COMPARISON WITH 
SEASAT ALTIMETER ANOMALIES 

* 4-LASER dataset is laser data from GEOS-1 , GEOS-2 , GEOS-3 and BE-C satellites 
3-LASER dataset is laser data from Dl-C , Dl-D , and PEOLE satellites 
OPTICAL dataset is the camera data from 20 satellites shown in TABLE 1A and IB 
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TABLE 4 

SUMMARY OF RESULTS FOR ERROR CALIBRATION 

CALIBRATION 

FACTOR 

GEM-T1 vs. GEM-T1 minus DATA SUBSET 

4-LASERS (GEOS 1,2,3 BE-C) 1.06 

STARLETTE LASER 1.10 

OSCAR + SEASAT DOPPLER 1.09 

OPTICAL ( 1 1 S ATS) 0.84 

LAGEOS LASER 1.45 

GEM-T1 vs. GEM-T1 + SURFACE GRAVITY 0.95 

GEM-T1 vs.GEM-Tl + SURFACE GRAVITY + 

SEASAT ALTIMETER 0.94 

GEM-T1 vs. SURFACE GRAVITY + SEASAT 

ALTIMETER 0.99 

GEM-T1 minus LAGEOS vs. LAGEOS 4- SURFACE 

GRAVITY + SEASAT ALTIMETER 0.95 

Weighting Factor f = 0.2 

10 TIMES DATA WEIGHT OF GEM-T1 

GEM-T1 vs. GEM-T1 minus 4-LASERS 2.75 
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Laser Residuals 



50 Points 

rms (r) * 7 a 


WEIGHT PER POINT « 1 ^ 

DU / \ <— 

(rms) 

.02 

“ i s2 

(rms ) 

.003 

“ 2 
a 

Figure 1. Characteristics of a Pass of Orbital Laser Residuals at a 
Tracking Site in Post Fit Analysis 
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RESIDUALS (MGAL 2 ) 


120 


100 H 


80 H 


60 H 


40 H 


30 H 


GEM-T1 WITHOUT 
KAULA CONSTRAINT 
(25X25) 



GEM-T1 


20 H 


* PGS3167 = GEM-T1 DATA SOLVED TO GEM-L2 SIZE 


I 

12 


16 


I 

20 


I 

24 


! 

28 


1 

32 


36 


DEGREE TRUNCATION IN HARMONICS 


Models show three modes of solution. The 25 x 25 field solves GEM-T1 tracking data without 
the Kaula constraint showing misclosure for high-degree terms. PGS-3167 solves GEM-T1 data 
(with Kaula constraint) to the GEM-L2 size field (20 x 20), showing no improvement over our 
previous model. GEM-T1 uses the Kaula constraint with a high-degree field (36x36) and is 
free of the above problems. 


Figure 2. Gravity Model Comparison With 1114 
5° x 5° SEAS AT Gravity Anomalies 
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